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Abstract. We present BVR photometric observations of the blazar PKS 0736+017. These observations were carried 
out with three telescopes in Mexico and two in Spain between December 1998 and April 2003. PKS 0736+017 
shows remarkable variation at different timescales and amplitudes. Maximum brightness was detected on December 
19, 2001 (B=14.90±0.01, V=14.34±0.01, and R=13. 79+0.01). A peculiar tendency to redden with increased 
brightness was detected throughout our observations. Moreover, in one season a good correlation between flux 
level and spectral slope is shown. This anomalous behaviour cannot be described by common flare models of 
blazars. The flux vs. spectral slope correlation observed in this and other blazars is worth further study. 
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1. Introduction 

Blazars, as a class, exhibit flux variability at timescales 
from minutes to years (e.g. Ghosh et al. 2000 Vagnetti 
et al. I20U51 Garcia et al. HMft Brown et al. 11989a! 
I1989bl Wright, McHardy & Abraham Ir39g|) . have high 
and variable polarization (Angel & Stockman I1980f) . and 
core dominated radio structures with weaker extended ra- 
dio emission (Antonucci & Ulvestad I1985fl . The contin- 
uum emission is characterized by a flat radio and steep 
optical-infrared spectrum, and the emission may extend 
to gamma-ray energies (Maraschi et al. 1994). Rclativistic 
beaming of the jet emission in radio-loud active galac- 
tic nuclei is generally accepted as the explanation for the 
observed properties of blazars (Blandford & Rees I1975|l . 
Their strong highly variable polarization and fluctuating 
brightness are usually taken to reflect processes at the 
heart of the nuclear engine. 

Some scenarios have been proposed to explain the 
rapid variability. Relativistic shock waves interacting with 
the inhomogeneous medium of the jet of blazars (Maraschi 
et al. HTM?! Qian et al. ITMT1 Marscher, Gear & Travis 
Marscher UnMI Gear et al. 1986; Ghisellini, Celotti 
& Costamante 2002) have been invoked; additional models 
refer to flares generated on the accretion disk of a super- 
massive black hole or to eclipsing of hot spots by the ac- 
cretion disk (Mangalam & Wiita 1993; Wiita lT§5oJl . The 
size of the emitting region is calculated using causality 
arguments. 



The amplitudes of variability have been found to be 
correlated with the redshift, and anticorrelated with lumi- 
nosities of AGNs (Hook et al. IHMl Cristiani et al. ITrMl 
Di Clemente et al. 119960 . With regard to timescales, no 
correlation was found with the amplitude of the variation 
in blazars by Nair l|1997f) or with luminosity or redshift in 
AGNs (Cristiani et al. I1996|l . An anticorrelation between 
timescales and luminosity of radio-loud quasars was found 
by Netzer et al. ifT^nTl . 

Many authors have found that the amplitude of the 
variations of blazars is higher at high frequencies, lead- 
ing to a flattened spectrum when the flux increases, and 
to a steeper one when it decreases (Racine 119701 Gear, 
Robson & Brown 119861 Massaro et al. 119981 Ghisellini 
et al. 119971 Maesano, Massaro & Nesci 119971 Maesano, 
Montagni, Massaro & Nesci ITW7|l . However, Ghosh et al. 
(2000) suggest that it may not be correct to make such 
generalizations. They found occasionally that the ampli- 
tudes of the variations are similar in V and R, in the ob- 
ject PKS 0235+164, or even larger in R, in the object 
PKS 0735+178. 

PKS 0736+017 was classified as a blazar by Angel 
& Stockman l|1980[l . and shows strong permitted and 
forbidden emission lines in the optical domain (Falomo, 
Scarpa & Bersanelli ©94). Padovani & Urry l|1992|l and 
Wall & Peacock (PSS) reported PKS 0736+017 as a 
Flat Spectrum Radio Quasar (FSRQ), while Brown et 
al. I|1989a[) report it as an Optically Violently Variable 
quasar, located at the centre of an elliptical galaxy, pos- 
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sibly a cD galaxy (Wright et al. 113551 Wyckoff, Gehren & 
WehingerUnEU Smith et al. 115551 Hutchings[T5J3 Falomo 
& Ulrich 125551 Hughes et al. 125551 Kotilainen, Falomo 
& Scarpa IT998[1 . which has two close resolved faint com- 
panions that are embedded in the nebulosity of the ob- 
ject (Dunlop et al. [X993; Huthchings, Johnson & Pyke 
115551 Bondi et al. H555)l . The absolute magnitude of the 
host galaxy has been determined between Mr=-24.3 and 
M R =-22.2 (Falomo & Ulrich 125551 Hutchings 1155?)) . and 
My=-21.7 (Smith et al. H555| . 

Historically, the brightness of this object has been 
variable at different timescales. Pica et al. 119801 I1988[) 
reported a 0.7 mag fade in six years and two flares of 
as much as 0.7 mag in 1974 and 1978, at timescales of 
weeks. Ghosh et al. (|2555|l observed PKS 0736+017 in 
the R band. They report variations of small amplitude 
(AR=0. 34+0.05) over timescales of six days. They also re- 
port a variation of 0.65+0.05 mag during the interval of its 
observations (between November 1995 and March 1998). 
Garcia et al. (1999) reported a variation of 0.8 mags (in V 
filter) in 400 days. McGimsey et al. Ijl975(l observed optical 
variation of 0.5 mag on a timescale of weeks and months, 
with an overall range of about 1.0 mag. Katajainen et al. 

(2000) detected brightness variations in V between 16.25 
and 16.84 magnitudes on timescales of 18 days. Schaefer 

(1980) reported Am B ^0.5 mags in 1 day. Ramirez et 
al. I)2003|l reported a variation of ~ 0.3 mag in 4.5 hours 
in the optical band. Recently, Clements, Jenks & Torres 
( 2003 ) found a flare with AR=1.3 mags in two hours early 
in 2002. 

Pica et al. lfT555)l detected that PKS 0736+017 was at 
its minimum in 1981, but immediately increased, reaching 
its peak brightness in 1984. These authors reported an 
average magnitude of 16.8 at the faint and 15.45 at the 
bright stages, with an average magnitude of 16.05 (in mp). 
The average magnitude for this object is between 16 and 
17 mag in the V band (e.g. Katajainen et al. 120001 Garcia 
et al. lT555l Wrig ht et a l. lT551 Bondi et al. lT55S|l . 

Bondi et al. (19921) detected variability at radio fre- 
quencies at two different epochs (in early 1980 and late 
1981). They found that the source was structured by 
two components and that the variability was caused by 
changes in both the flux densities and the separation of the 
components. In a third epoch, they identified the bright- 
est component with the core. McAdam fI976$ reported 
flux variability at 408 MHz, and Kedziora-Chudczer et al. 

( 2001 ) observed polarization variations. 

On the other hand, other studies were carried out with- 
out positive results. Padrielli et al. (1987) made multifre- 
quency observations of PKS 0736+017 to search for cor- 
relation between low and high frequencies. Jones et al. 

( 1981) observed it as a part of their study of twenty radio 
sources to try to categorize flux variations due to changes 
in source scale, structure or the slope of the electron en- 
ergy spectrum. Sitko, Schmidt & Stein (|1985f) searched for 
optical polarimetry. 

Brown et al. (|1989b|) . within their investigation of the 
variability properties from centimetre to ultraviolet wave- 



lengths, reported that they did not find any correlation 
between near-infrared flux level and spectral slope for 
PKS 0736+017. However, they noted that excluding the 
data for the minimum, the flux shows a significant anti- 
correlation (more than 80% confidence level) between flux 
level and spectral slope (taken /„ cx v a ). 

In this paper we present and discuss optical vari- 
ability of the Optical Violently Variable (OVV) quasar 
PKS 0736+017; in §2 the observations and data reduction 
procedures are described; in §3 the variability behaviour 
is analyzed. Finally, in §4 we present our discussion and 
in §5 our conclusions. 

2. Observations and data reduction 

The observations were carried out with five telescopes in 
Mexico and Spain. Observations from Mexico were made 
with three telescopes: a 2.1 m (Ml), a 1.5 m (M2), and a 
0.84 m (M3) operated by the Observatorio Astronomico 
Nacional at San Pedro Martir, Baja California. Telescopes 
in Spain were the 1 m Jacobus Kapteyn Telescope (SI) op- 
erated by the Isaac Newton Group at the Observatorio 
del Roque de los Muchachos, and the 1.52 m (S2) 
Estacion de Observation de Calar Alto (EOCA) at Centro 
Astronomico Hispano-Aleman Calar Alto (CAHA) oper- 
ated by the Observatorio Astronomico Nacional. All the 
observations were made with BVR Johnson series filters. 
The log of the observations is given in Table 1 . In this table 
the date and telescope used are listed in columns (1) and 
(2), respectively. Magnitudes and uncertainties for each 
filter are listed in columns (3)- (8). 

The CCDs used with Ml and M2 were a Site 
SI003 with 1024X1024 pixels of 24/xto x 24/im with 
a Metachrome II and VISAR coating to improve the 
response in the blue. With M3 the detector was a 
Thomson TH7398M CCD with 2048X2048 pixels, of 
14/irn x 14/im, with a Metachrome II coating. SI had at- 
tached a TEK1024AR: 1024X1024 pixels of 24/xto x 24/ito, 
thinned AR coated, and S2 had a Tektronics TK1024AB: 
1024X1024 pixels of 24/xto x 24/iTO. 

The results of five observing seasons (one was carried 
out with two telescopes, M2 and S2) are reported here. 
Four nights out of fifteen were rejected due to observing 
problems like extremely bad seeing and clouds. From the 
remaining eleven nights of observation, one was taken with 
the Ml, two with the M2, five with the M3, one with 
the SI, and two with the S2. On these data we detected 
year-to-year, month-to-month, day-to-day, and intra-night 
variability. Only non- variable stars in the same frame were 
taken into consideration. 

The observations were done with the strategy ex- 
plained by de Diego et al. (1998), taking five images with 
each filter in BVR sequence with approximately one hour 
between each observation set. The observations were re- 
duced with the IRAF-APPHOT package. The best aper- 
ture was determined using the growth curves technique, 
which resulted in a pseudoaperture ~3 FWHM. Standard 
stars given by Landolt (1953) were used to obtain our mag- 
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Table 1. Summary of observations. See the text for a 
full description of the table. The magnitudes listed for 
December 2001 are a mean values. 



Date 


Tels. 


m,B 


Stub 


mv 


Smv 


VTbR 


SrriR 


Dec/11/1998 


SI 


16.91 


0.02 


16.50 


0.01 


15.97 


0.02 


Afar/03/2000 


M2 


16.60 


0.04 


16.09 


0.03 


15.56 


0.03 


Mar /04/2000 


S2 


16.61 


0.01 


16.10 


0.01 


15.60 


0.01 


Mar /05/2000 


52 


16.29 


0.02 


15.78 


0.01 


15.25 


0.03 


7Vou/26/2000 


M3 


17.20 


0.02 


16.70 


0.01 


16.17 


0.01 


7Vcw/28/2000 


M3 


17.29 


0.02 


16.83 


0.02 


16.32 


0.02 


TVcw/29/2000 


M3 


17.30 


0.02 


16.80 


0.02 


16.26 


0.02 


Dec/03/2000 


M3 


17.19 


0.02 


16.63 


0.03 


16.08 


0.02 


Dec/04/2000 


M3 


16.98 


0.02 


16.48 


0.02 


15.89 


0.02 


Dec/19/2001 


Ml 


15.07 


0.01 


14.49 


0.01 


13.94 


0.01 


Apr /26/2003 


M2 


16.90 


0.03 


16.36 


0.03 


15.74 


0.03 



nitudcs. Uncertainties were determined from magnitude 
dispersion within each image sequence, which docs not 
differ significantly from dispersion between the groups of 
images for the reference stars. When the data were trans- 
formed to fluxes, corrections to the rest frame of the object 
and dereddening were carried out. 

We detected three timescales of variability: one 
year (~2.3 mag in BVR, between November 2000 and 
December 2001), one day (^0.1 in BVR, within the 
November 2000 season), and intra-night (~0.3 mag in 
BVR, in December 2001). These variations are studied 
in the next section. 



3. Analysis of variability 

3.1. Variability 

PKS 0736+017 was observed within a larger program to 
study AGN microvariability. A few observations were also 
made to monitor long-term variability. Thus, the data per- 
mit us to obtain relevant information about short- and 
long-term variability. However, the sparse long-term sam- 
pling limits the possibilities for determining parameters 
such as the precise moment when an outburst started and 
finished, or the date on which the maximum occurred. 

Figures 1 and 2 show the BVR light curves, which show 
flux variations with timescales of years and minutes, re- 
spectively. The long-term variations have amplitudes up 
to a few magnitudes, while the amplitudes for the short- 
term variability range from a few hundredths to several 
tenths of a magnitude. In spite of the poor time sampling, 
the continuous activity is obvious. Fig. 1 also shows light 
curves for three of our control field stars. Some stars in 
the PKS 0736+017 field show variability of much smaller 
amplitude (less than 0.03 mags), but only those that re- 
main constant during the observations were selected. In 
addition, star 2 from Smith et al. ©85) for the field 
of PKS 0736+017 was observed. Our magnitude estima- 
tions for this star (B=15.92, V=15.40, and R=15.04) dif- 
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Fig. 1. The left panel shows the light curves in B,V, and 
R bands for PKS 0736+017, as listed in Table 1. The right 
panel shows the reference stars included in our field: down 
triangles for star 9, squares for star 10, and up triangles 
for star 11. The December 2001 points are a mean value. 
J.D. = 2451158.5 corresponds to December 11, 1998. 



fer slightly from those given by Smith et al. (B=15.85, 
V=15.37, and R=14.97). 

An important event was detected on December 19, 
2001, when the flux level reached a maximum of B~14.90, 
V~14.34, R^13.79 (in Tables 1 and 2, and Fig. 1 we 
only show mean values of the magnitudes for this night). 
Because of the sparse sampling, the real maximum is 
not well timed, but it should have happened around this 
epoch (see Fig. 1). Added to this outburst, PKS 0736+017 
showed a microvariability event during the 4.5 hours that 
it was monitored (Fig. 2c). It decreased its brightness at 
a rate of ^0.066 magnitudes/hour. Due to instrumental 
problems we lost the first R point. 

Previous photometric observations showed that this 
blazar varied between 16 and 17 mag in V (e.g. Garcia et 
al.EESI Bondi et aL lUMI Wright et al.HSSHI Katajainen 
et al. I2000[) . Throughout our observations the brightness 
was generally between these historical values. In December 
2001, we detected the brightest magnitude (see Table 1). 

In general, the amplitudes of the variations are much 
larger than the uncertainties, and statistical techniques 
are unnecessary (the smallest detected flux variations, be- 
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Fig. 2. (a-c).a)Light curves for the March 2000 season, 

and b) for November 2000. Magnitudes for December 2001 TaWe 2 Comparison of the spectral slope relative to the 

shown in c are instrumental. n ., , 

K magnitude. 



tween November 26 and 28, 2000, are 0.09±0.03 magni- 
tudes in B, 0.13±0.02 magnitudes in V, and 0.15±0.02 
magnitudes in R). 

We presume that this activity began at least in early 
1998, since Katajainen et al. (2000) reported V magni- 
tudes in the range from 16.84 to 16.25 between 1995 and 
1997. Later, Ghosh et al. (2000) observed an R magnitude 
of 16.24 (within the typical values) on December 15, 1997, 
and on March 25, 1998 they detected that the brightness 
had increased to 15.47, and finally diminished to 15.94 in 
April 1998 in the same band. On December 11, 1998 we 
detected that the brightness had dimmed to V=16.50, and 
R=15.97. 

The object was monitored to study its long-term vari- 
ability in December 1998 and April 2003. Its brightness 
was similar in both epochs (V ~ 16.4), although the ob- 
served colour was different, bluer in December 1998 than 
in April 2003 (see section 3.2 and Table 2). 

3.2. Variability behaviour 

The most remarkable behaviour through all our obser- 
vations is that PKS 0736+017 became redder when it 
was brighter (see Figs. 3, 4 and 5). In other words, the 
spectrum became steeper when the object was bright and 
flatter when it faded. This behaviour is opposite to the 
general reports for this class of objects (e.g. Kedziora- 
Chudczer et al. l2lM| Gear et al-HD Racine 1970; Gear, 



Date ran ct 5a 

Dec/11/1998 15.97 1.24 0.05 

Mar/03/2000 15.56 1.40 0.07 

Mar/04/2000 15.60 1.38 0.02 

Mar/05/2000 15.25 1.45 0.04 

iVo«/26/2000 16.17 1.40 0.04 

iVo«/28/2000 16.32 1.30 0.04 

iVou/29/2000 16.26 1.41 0.05 

Dec/03/2000 16.08 1.58 0.07 

Dec/04/2000 15.89 1.55 0.05 

Dec/19/2001 13.94 1.72 0.02 

Apr/26/2003 15.74 1.67 0.02 



Robson & Brown IT§g6l Massaro et al. Maesano et 

al. Q997 Papadakis et al. I2()()3[) . and common models 
(e.g. Li and Kusunose IMtffl Chiang and Bottcher IMHSl 
Wang and KusunoseEEEl Dermer lUMl Spada et al. lMfTl 
Sikora et al lSUUTl Qian et al. lTMTl Wiita HUMI Dermer & 
Schlickeiser l2Hrj2l Nesci et al- HMSI Marscher et al. H^SUI 
Trevese & Vagnetti |?002 ) which state that blazars be- 
come bluer when they brighten. This peculiar behaviour 
was also noted by Clements et al. (2O03J in early 2002, 
just a month after the flare of December 2001. These au- 
thors reported V=14.322 and R=13.717 (V-R ~0.605) at 
the maximum. They concluded that the phenomenon "ap- 
pears to be more related to the nature of the variation than 
to the host galaxy or the brightness level" . 
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Fig. 3. The shape of the optical broad band of PKS 
0736+017 can be fitted by a power law. 



1.50 



S 1.40 



1.30 




1.30 



1.25 



1.20 

0.10 0.20 0.30 -0.20 -0.13 -0.07 0.00 0.07 



log(f)(mJy) 



log(f)(mJy) 



Fig. 4. a) Data for March 2000 shows a tendency in the 
sense that PKS 0736+017 is redder when it is brighter, 
b) This tendency was confirmed with data obtained in 
November 2000, which shows a clear correlation between 
colour and brightness. 



Our data can be described well with a power law, i.e., 
f„ oc v~ a , where —a is the spectral slope in a log flux vs. 
log frequency diagram (see Fig. 3), which permits us to 
quantify the reddening accurately. Better fits can be made 
if more spectral components are included in the spectrum 
model, such as the contribution of the galaxy (e.g. Wright 
et al-HSHJ or a blackbody (e.g. Malkan and Moore 11555)1 . 
but with spectral resolution of only three points an extra 
component does not make sense. 

The values of the spectral slope and their uncertainties 
for each epoch are given in columns (3) and (4) in Table 
2. In this table the date and the R magnitude are listed 
in columns (1) and (2), respectively. It is evident that 
commonly the brightest points have steeper spectra. This 
tendency was noticed in November 2000, when we detected 
a slope steepening from a =1.30 to 1.55 while the flux 
increased (Fig. 4b). Furthermore, the flux and the spectral 
slope are well correlated, with a correlation coefficient of 
89 %. In March 2000 we have two estimates of a (1.38 and 
1.45) and one of them with a large uncertainty. However, 
Fig. 4a shows that the tendency is still the same as in Fig. 
4b. At this time, the brightness increased while the slope 
changed slightly. Note that the uncertainties for March 
2000 do not overlap. 



A similar tendency to redden with increased bright- 
ness was observed in the long-term variability (see Fig. 
5) , but a marginal correlation between spectral index and 
flux level was found. Some data differ from the common 
behaviour in our observations. In particular, in December 
2001, while the flux decreased by 0.3 mags the spectral 
slope remained constant (a=1.72). On the other hand, 
April 2003 data has a low flux level and high spectral slope 
value. Nevertheless, we can find an explanation in the on- 
line data of Clements et al. ( 2003). On an a- log flux plane, 
these data are distributed in two varying modes: one pro- 
ducing larger flux variations with shorter spectral changes 
(tendency represented by dashed line in Fig. 5), and other 
one with the inverse effect (tendency represented by dot- 
ted line in Fig. 5). Then, December 2001 could be under- 
stood if it is varying as the last variable mode, while April 
2003 as the former one. 

In December 2001 both flux level and behaviour were 
similar to those reported by Clements et al. I|2003|l for 
the January 2002 flare. In Fig. 5, Clements et al.'s average 
maximum from January 23 is plotted (roughly computed 
with only R and V bands). 
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4. Discussion 

We stated that the spectrum can be described well by 
a power law, although other spectral components could 
be included. For example, Malkan and Moore (.1986) ob- 
served PKS 0736+017 at ultraviolet, optical and infrared 
bands. In the optical and ultraviolet bands they fitted a 
power law with spectral index of 1.0 plus a blackbody 
with a temperature of 26,000 K, but the latter should only 
contribute a fourth of the total brightness. Additionally, 
these authors suggested that PKS 0736+017 has combined 
emission from a "normal" quasar (with a mixture of both 
power law emission with a spectral slope of 1.1, and ther- 
mal emission from accretion disk) and a blazar emission 
(which is represented by a power law with a slope between 
1.1 and 2.0). They further proposed that the blazar com- 
ponent should be increasingly important at longer wave- 
lengths, and when the continuum is brighter. 

Brown et al. (1989a) also proposed that the emission 
of OVV quasars is a mixture of a "quiescent" , a violently 
variable and a thermal emission component. The first com- 
ponent should be responsible for the centimetre emission, 
the second one dominates the millimetre to ultraviolet re- 
gion, and the third is responsible for some of the opti- 
cal and ultraviolet emission. With this picture in mind, 
the violently variable component should cause the flares 
and outbursts in these objects. Brown et al. ( 1989a 1989b) 
also mentioned that, at high brightness state, the spectral 
shape of OVV quasars showed similar properties to the 
emission of BL Lac objects. 

These simple models could give a qualitative explana- 
tion only for our low flux level data. This can be noticed 
considering that, the composite spectrum is flatter than 
the non-thermal component, because the thermal contri- 
bution is larger in the blue part. Then, when the object is 
brightening, the non-thermal component has an even more 
dominant contribution to the total flux, and the composite 
spectrum steepens. 

However, high flux levels have no satisfactory explana- 
tion because our point of largest slope, in December 2001, 
corresponds to an outburst (at this time, the blazar com- 
ponent must be dominant and the flattest spectra must be 
measured, as is observed in these objects). The bulk of our 
data differ from synchrotron loss energy behaviour (Brown 
et al. H989a|) . A thermal varying bump cannot explain our 
data because changes of this component lead to larger 
changes in the blue part of the spectrum. Additionally, 
two variability modes in the a-log flux plane are observed. 

We have detected a behaviour that is not explained by 
the usual mechanisms. The evidence that the amplitudes 
of variations are not systematically larger at high frequen- 
cies in this and other blazars is increasing (see for example 
Ghosh et al.EOOO Brown et al. H989bl Malkan and Moore 
HM61 Massaro et al. ITMSl Clements et al. IMiHl and this 
work) . 

Brown et al. I|1989b|> observed PKS 0736+017 at in- 
frared wavelengths and found no correlation between spec- 
tral slope and flux level; however, if they disregard one 
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Fig. 5. Relation between the spectral index a and the 
brightness in R. We have included our own data and 
the average maximum reported by Clements et al. (2003) 
for January 2002. Our analysis of the data reported by 
Clements et al., in their article online version, suggests 
two varying modes. This tendency is represented by the 
dashed and dotted lines. Our data agrees with these ten- 
dency. At low flux level, small changes in flux correspond 
to large changes in the spectral slope, while much less pro- 
nounced spectral changes correspond to high brightness 
state. In both cases, the object reddens when it brightens. 



point at the minimum, they find that the spectrum steep- 
ens when it is brighter, as they also found in the case of the 
blazar 3C 279. They did not find a significant correlation 
between the flux level and the slope values for other three 
blazars (1308+326, 3C 345, and 3C 446). They point out 
that this absence of correlation is not consistent with the 
behaviour expected from emission of a single synchrotron 
component, therefore other components must be involved. 

Earlier, Miller lfTU5T|l found that the blazar 3C 446 be- 
came redder as the source brightened. Brown et al. Q1986[) 
reported that the spectrum of 3C 446 steepened with 
brightening in the K-band in July 1980 and September 
1983. They explained their observations with the argu- 
ment that the maximum of brightness could have occurred 
before their observations and that the spectrum steepened 
rapidly. This explanation cannot be applied to our data 
because we observed PKS 0736+017 when it was bright- 
ening (November 2000). 
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Ghosh ct al. ( 2000 ) found a reddening in their optical 
observations of the BL Lac object PKS 0735+178. 
Additionally, they observed that in the blazar 
PKS 0235+164, the spectral slope remained almost 
constant when its brightness increased. They concluded 
that these characteristics cannot be described by energy 
losses in a pure synchrotron mechanism scenario (one of 
the most used flare models). 

5. Conclusions 

The OVV quasar PKS 0736+017 was observed between 
December 1998 and April 2003. Variability events were de- 
tected throughout the observations at different timescales 
and amplitudes. Some flares were detected in March 2000 
and November-December 2000, with flux change of ^0.32 
and ~0.4 mags, respectively. 

The most important features that we report are: 

• An outburst was detected on December 19, 2001. At 
this time, the bright level reached a maximum of 14.90 in 
the B band, 14.34 in the V band, and 13.79 in the R band. 
A brightness decrease of ~0.3 mags occurred in 4.5 hours, 
while the spectral slope was unchanged. 

• A tendency to redden with increased brightness was 
detected throughout our observations. Moreover, a good 
correlation between flux level and spectral slope was de- 
tected in November 2000, while a clear tendency was 
found in March 2000. Clements et al. (2003) present data 
strongly suggesting the existence of two varying modes. 
Our data agree with this assumption. 

Clements et al. ()20031 also detected the same tendency 
between brightness and spectral slope for PKS 0736+017, 
but other authors did not (e.g. Brown et al. H989bl Malkan 
& Moore I1986f) . There is no doubt that more observa- 
tions of blazars that show reddening as they brighten 
are needed. Multiwavelength observations with good time 
resolution on targets like 3C 446, PKS 0735+178, and 
PKS 0736+017 will give us important information about 
the physical conditions of these objects. Models should be 
improved to account for these observational results. 

Acknowledgements. This work was supported by Grant No. 
121553 from CONACyT, and PAPIIT projects IN115599, and 
IN113002. The Jacobus Kapteyn Telescope is operated on the 
island of La Palma by the Isaac Newton Group in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de 
Astrofi'sica de Canarias. We acknowledge the use of IRAF, dis- 
tributed by NOAO. We also acknowledge to Garcia Ruiz, G, 
Melgoza, G, Monrroy, S., and Montalvo, F. for their support 
during the observations. 

References 

Angel, J.R.P. & Stockman, H.S. 1980, ARA&A, 18, 321 
Antonucci, R.R.J. & Ulvestad, J.S. 1985, ApJ, 294, 158 
Blanford, R.D. & Rees, M.J. 1978, in Wolfe A.M. ed. 

Pittsburgh Conference on BL Lac Objects. University of 

Pittsburgh, p. 328 
Bondi, M., Padrielli, L., Fanti, R. et al. 1996, A&A, 308, 415 



Brown, L.M.J. , Robson, E.I., Gear, W.K. et al. 1986, MNRAS, 
219, 671 

Brown, L.M.L., Robson, E.I, Gear, W.K. et al. 1989a, ApJ, 
340, 129 

Brown, L.M., Robson E.I., Gear, W.K. & Smith, M.G 1989b, 

ApJ, 340, 150 
Chiang, J. & Bottcher, M. 2002, ApJ, 564, 92 
Cristiani, S., Trentini, S., La Franca, S. et al. 1996, A&A, 306, 

395 

Clements, S.D., Jenks, A. & Torres, Y. 2003, AJ, 126, 37 

de Diego, J. A., Dultzin-Hacyan, D., Ramirez, A. & Bem'tez, E. 

1998, ApJ, 500, 69 
Dermer, CD. 1998, ApJ, 501, L157 
Dermer, CD. & Schlickeiser, R. 2002, ApJ, 575, 667 
Di Clemente, A., Giallongo, E., Natali, G, Trevese, D. & 

Vagnetti, F. 1996, ApJ, 463, 466 
Dunlop, J.S., Taylor, G.L., Hughes, D.H. & Robson, E.I. 1993, 

MNRAS, 264, 455 
Falomo, R., Scarpa, R. & Bersanelli, M. 1994, ApJS, 93, 125 
Falomo, R. & Ulrich, M.H. 2000, A&A, 357, 91 
Garcia, A., Sodre, L., Jablonski, F.J. & Terlevich, R.J. 1999, 

MNRAS, 309, 803 
Gear, W.K., Brown, L.M.J. , Robson, E.I. et al. 1986, ApJ, 304, 

295 

Gear, W.K., Robson, E.I. & Brown, L.M.J. 1986, Nature, 324, 
546 

Ghisellini, G, Celotti, A., Costamante, L. 2002, A&A, 386, 833 
Ghisellini, G, Villata, M., Raiteri, CM. et al. 1997, A&A, 327, 
61 

Ghosh, K.K., Ramsey, B.D., Sadun, A.C. & Soundararaja- 

perumal, S. 2000, ApJS, 127, 11 
Hook, I.M., McMahon. R.G, Boyle, B.J. & Irwin, M.J. 1994, 

MNRAS, 268, 305 
Hughes, D., Kukula, M.J., Dunlop, J.S. & Boroson, T. 2000, 

MNRAS, 316, 204 
Hutchings, J.B. 1987, ApJ, 320, 122 

Hutchings, J.B., Johnson, I. & Pyke, R. 1988, ApJS, 66, 361 
Jones, W.T., Rudnick, L., Owen, F.N. et al. 1981, ApJ, 243, 
97 

Katajainen, S., Takalo, L.O., Sillampaa, A. et al. 2000, A&A, 
143, 357 

Kedziora-Ctmdczer, L.L., Jauncey, D.L., Wieringa, M.H., 

Tzioumis, A.K. & Reynolds, J.E. 2001, MNRAS, 325, 1411 
Kotilainen, J., Falomo, R. & Scarpa, R. 1998, A&A, 332, 503 
Landolt, A.U. 1992, AJ, 104, 340 
Li, H. & Kusunose, M. 2000, ApJ, 536, 729 
Maesano, M., Messaro, E. & Nesci, R. 1997, IAU, Circ. 6700 
Maesano, M., Montagni, F., Massaro, E. & Nesci, R. 1997, 

A&A, 122, 267 
Malkan, M. & Moore, R. 1986, ApJ, 300, 216 
Mangalam, A.V. & Wiita, P.J. 1993, ApJ, 406, 420 
Maraschi, L., Grandi, P., Urry, CM. et al. 1994, ApJ, 435, L91 
Maraschi, L., Celotti, A. & Treves, A. 1989, in Proc. 23d, 

ESLAB Symp. on Two Topics in X-ray Astronomy ed. J.J. 

Hunt & B. Battrick (vol. 2; Paris: ESA), 825 
Marscher, A.P. 1996, in ASP Conf. Ser. 110, Blazar continuum 

Variability, ed. H.R. Miller, J.R. Webb & J.C Noble (San 

Francisco: ASP), 248 
Marscher, A. P., Gear, W.K. & Travis, J.P. 1992, in Variability 

of Blazars, ed. E. Valtaoja & M. Valtonen (Cambridge 

Univ. Press), 85 
Marscher, A.P. 1980, ApJ, 235, 386 

Massaro, E., Nesci, R., Maesano, M., Montagni, F. & 
D'Allessio, F. 1998, MNRAS, 299, 47 



8 



Ramirez et al.: Optical variability of PKS 0736+017 



McAdam, W.B. 1976, PASA, 3, 86 

McGimsey, B.Q., Smith, A.G., Scott, R.L., et al. 1975, AJ, 80, 
895 

Miller, H.R. 1981, ApJ, 244, 426 
Nair, A.D. 1997, MNRAS, 287, 641 

Nesci, R., Maesano, M., Massaro, E. et al. 1998, A&A, 332, LI 
Netzer, H., Heller, A., Loinger, F. et al. 1996, MNRAS, 279, 
429 

Padovani, P. & Urry, CM. 1992, ApJ, 387, 449 

Padrielli, L., Aller, M.F., Aller, H.D. et al. 1987, A&AS, 67, 63 

Papadakis, I.E., Boumis, P., Samaritakis, V. & Papamastora- 

kis, J. 2003, A&A, 397, 565 
Pica, A. J., Pollock, J.T., Smith, A.G. et al. 1980, AJ, 85, 1442 
Pica, A.J., Smith, A.G., Webb, J.R. et al. 1988, AJ, 96, 1215 
Qian, S.J., Quirrembach, A., Witzel, A. et al. 1991, A&A, 241, 

15 

Racine, R. 1970, ApJ, 159, L99 

Ramirez, A., de Diego, J. A., Dultzin-Hacyan, D. & Bemtez, 
E., 2003, in "Active Galaxy Nuclei: From central engine to 
Host Galaxy", ASP Conference Series, Vol. 290, 2003, eds. 
S.Collin, F. Combes, and I. Shlosman, p. 123 

Schaefer, B.E. 1980, PASP, 92, 255 

Sikora, M., Blazejowski, M., Begelman, M.C. & Moderski, R. 

2001, ApJ, 554, 1 
Sitko, M.J., Schmidt, G.D. & Stein, W.A. 1985, ApJS, 59, 323 
Smith, E.P., Heckman, T.M., Bothun, G.D., Romanishin, W. 

& Balick, B. 1986, ApJ, 306, 64 
Smith, PS., Balonek, T.J., Heckert, PA. et al. 1985, AJ, 90, 

1184 

Spada, M., Ghisellini, G., Lazzati, D. & Celotti, A. 2001, 

MNRAS, 325, 1559 
Trevese, D. & Vagnetti F. 2002, ApJ, 564, 624 
Vagnetti, F., Trevese, D., & Nesci, R. 2003, ApJ, 590, 123 
Wall, J.V. & Peacock, J.A. 1985, MNRAS, 216, 173 
Wang, J.-M. & Kusunose, M. 2002, ApJS, 138, 249 
Wiita, P.J. 1996, in ASP Conf. Ser. 110, Blazar continuum 

Variability, ed. H.R. Miller, J.R. Webb & J.C. Noble (San 

Francisco: ASP), 42 
Wright, A.E. 1984, PASA, 5, 510 

Wright, S.C., McHardy, I.M. & Abraham, R.G. 1998, MNRAS, 
295, 799 

Wyckoff, S., Gehren, T. & Wehinger, P.A. 1981, ApJ, 247, 750 



